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We calculate the non-equilibrium charge transport properties of nanoscale junctions in the steady 
state and extend the concept of charge susceptibility to the non-equilibrium conditions. We show 
that the non-equilibrium charge susceptibility is related to the non-linear dynamical conductance. 
In spectroscopic terms, both contain the same features versus applied bias when charge fluctuation 
occurs in the corresponding electronic resonances. However, we show that while the conductance 
exhibits features at biases corresponding to inelastic scattering with no charge fluctuations, the 
non-equilibrium charge susceptibility does not. We suggest that measuring both the non-equilibrium 
conductance and charge susceptibility in the same experiment will permit us to differentiate between 
different scattering processes in quantum transport. 
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Recent developments in modern techniques of micro- 
scopic manipulation and nanotechnologies enable us to 
build functional nanoscale systems, for example, elec- 
tronic nanodevices or molecular motors llitSj- In such 
systems, classical equilibrium thermodynamics is not 
suited to describe the quantization of the charge or heat 
flow. The properties of such devices differ from their 
equilibrium counter parts since non-equilibrium quantum 
and non-linear effects dominate. The concepts of con- 
ventional statistical mechanics and linear-response the- 
ory for these small systems need to be substituted with 
those of non-equilibrium quantum statistical mechanics 
[^ftSj ■ This is the appropriate framework for dealing with 
nanoscale systems driven out of equilibrium, especially 
when one wants to design or control these systems as 
heat engines or electro-mechanical devices. 

Extending the concept of equilibrium statistical me- 
chanics (linear-response theory and response functions, 
fluctuation theorems) to the non-equilibrium (NE) con- 
ditions has seen a recent growing interest It is 
known from linear-response theory that there exists some 
relationship between different response functions, like for 
example, the density-density, current-density, or current- 
current correlation functions at equilibrium 'TT . How- 
ever there is no reason why these relationships should 
hold at NE. Motivated by understanding these NE prop- 
erties and their use in practical nanoscale devices, we 
focus in this paper on a specific physical property: the 
electronic transport. In particular, we consider the re- 
lationship between the electrical conductance and the 
charge susceptibility in nanoscale junctions. We provide 
a definition for the NE charge susceptibility, which can 
be measured in experiments, and examine in detail its 
relationship with the full non-linear dynamical conduc- 
tance. 

We show that the non-equilibrium charge susceptibil- 
ity and the dynamical conductance of such a system are 



related to each other, though in a different manner than 
at equilibrium. At finite bias, they both contain informa- 
tion about the charge fluctuation (induced by the bias) 
in the electronic resonances. However the NE charge sus- 
ceptibility does not contain information about purely in- 
elastic scattering processes which do not involve charge 
fluctuations. By measuring both the conductance and 
the NE charge susceptibility in the same experiment, one 
can identify the nature of scattering processes involved 
in transport through nanoscale junctions. 

We illustrate this property with numerical calcula- 
tions for a model of a single-molecule nanojunction 
in the presence of electron-phonon coupling. Our re- 
sults are relevant, but not limited only, for electron- 
phonon scattering processes. Other examples could be 
electron-plasmon, electron-electron, electron-spin excita- 
tions scattering events. In the following, we first briefly 
recall the relationship between linear conductance and 
charge susceptibility at equilibrium. Then we derive 
the corresponding relationship in the NE conditions, and 
present numerical calculations. 

Equilibrium response functions: Within the 
linear-response theory of a system at equilibrium |12[ 1141 
I15j . the current / is related to a frequency-dependent 
applied bias V via the linear conductance g as I{oj) = 
g{uj)V{uj). The linear conductance is a response func- 
tion obtained from the current-density correlation func- 
tion g{t) = {ie/h){[i{t),N{Q)])e{t), where N is the to- 
tal occupancy operator and / is the current operator 
/ — edN /At. The linear conductance g is directly re- 
lated to the density-density correlation function Xcit) = 
-i{[N{t),N{Q)])6{t) by the relation g{uj) = iw^Xcl^). 
Xc is also known as the charge susceptibility and repre- 
sents the response function of the charge density modi- 
fications 5n due to variation of the electrostatic poten- 
tial 5v: Sn — J Xc5v [16]. In the DC limit, one gets a 
finite linear conductance g{uj — > 0) when the charge sus- 
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ceptibility goes as Xc{^) = f{^)/^ '^ith /(O) 7^ 0. At 
equilibrium, there is a clear and well-defined relationship 
between the charge susceptibility and the linear conduc- 
tance. However there is no a priori reason why such a 
relation should still hold at non-equilibrium when an ap- 
plied bias drives the system in a non-linear regime. 

Non-equilibrium charge susceptibility and 
transport: We consider a generic system consisting of 
a interacting central region C, the scatterer of inter- 
est (e.g. a molecule or a quantum dot), connected to 
two electrodes, acting as source and drain. The elec- 
trodes are non-interacting Fermi seas at their own equi- 
librium and there is no direct contact between them. We 
use non-equilibrium Green's functions (NEGF) to cal- 
culate the electric current and charge of the system in 
NE conditions [13]. The system is under a finite, but 
not small, applied bias and is assumed to have reached 
a non-equilibrium steady-state which can be described 
by an effective (pseudo) equilibrium steady-state density 
matrix [SjlTT]. 

We define the non-equilibrium charge susceptibility 
in the steady-state as the response (not necessarily 



ductance G{V) can be written as: 



x; 
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linear) for the modifications of the total electronic occu- 
pancy of the central region 5{nc) due to the changes in 
the applied bias 5V , i.e. changes in the cause that drives 
the system out-of-equilibrium [18] : 



dV 



(1) 



The total occupancy 



,NE\ 



of the central region C is 



given by the non-equilibrium lesser Green's function as 
(n^^) = — i/ da;Tr[G'^(aj)]/27r, where the trace runs over 
the electronic states in the region C. 

We now examine in detail how x^^(V^) is related to 
the dynamical conductance G{V) = dl /dV . The current 
at the left L interface between the central region C and 
the L lead is given by the Meir-Wingreen expression [13] : 



II = 



g Tr [/L(c.)(G^(a.) - G'h{u:))V l{u^) 
+G<{lo)Tl{u)\ , 



with Ti^{lo)/2 being the imaginary part of the L lead 
self-energy, and G^"' being the retarded, advanced and 
lesser Green's function of the central region respectively, 
and the trace is taken over the electron states of the cen- 
tral region. 

By using the properties of a NE steady state, one intro- 
duces a non-equilibrium distribution functional /^^ for 



the central region as GS,{lo) 



-f^^{io){G^c-G'h){^) 



[19j . At equilibrium is simply given by the conven- 
tional Fermi distribution function. The dynamical con- 



G{V) 



J^Tr[{l-fUfn-')dvG<T, 

-5y(/i(/n-i)G<r^], 



(3) 



which shows a relation between the dynamical conduc- 
tance and the derivative of the lesser Green's function 
versus the applied bias dyG^. To show more clearly how 
G{V) and x?^(^) are related to each other, we consider 
the following simpler system. 

A model system: The model consists of a single elec- 
tron level in the region G, in the presence of some arbri- 
tary kind of interaction. For the moment, we consider 
the wideband limit where Ti^(uj) = F, and that all the 
potential drop occurs at the left contact. Only the Fermi 
distribution Jl of the left lead depends explicitly on the 
bias V via the Fermi level ^l- Within these conditions, 
we find a relation between the dynamical conductance G 
and the non-equilibrium charge susceptibility x^^- 

G{V){^T)-' + xr(^^) = / dujdvihAci^)), (4) 

where Ac{uj) = (G^(a;) - Gj;(w))/27ri. 

For non-interacting systems, the spectral function Ac 
is independent of the bias, then dvAc{i-o) = 0. By using 
the definitions of G and for symmetric contacts and 
the corresponding non-equilibrium distribution function 
/^^ = {Tl/l + TM/iTL + Ffl) = iU + /fl)/2 m. 
we find a direct proportionality between G and x?^- 
G{V) = ^TXc^iV)/e [21]. Beyond the wideband ap- 
proximation (with symmetric contacts), we obtain the 
relation: G(F) = '^r{fiL)x^^{V)/e. Hence the com- 
patibility between the equilibrium and NE approaches 
implies that lim^_j.o iwXc('^) = rx?^(^) (within the DC 
limit of linear- response) . 

For interacting systems, Ac depends on V through 
the interaction self-energy Ei„t(w,y). An analytical ex- 
pression relating G and x?^ is more difficult to obtain 
[22] . However we show next, from numerical calculations 
beyond the wideband limit, that there is a clear relation- 
ship between G{V) and X?^(^) for a model of interaction 
self-energy. 

An application: For this we have to make a choice 
for the interactions in the central region C. The NE 
charge susceptibility has been briefly studied for a model 
of electron-electron interaction in the Anderson impu- 
rity model at non-equilibrium in [23j . In the following 
we consider a model electron-phonon interaction in the 
central region G [231 [5S] . Considering such a model per- 
mits us to get several different physical effects: the renor- 
malization of the electron level but also all the phonon 
replica (the phonon side-band peaks). So effectively we 
are dealing with a richer model of multi-electronic reso- 
nances. Such a model includes different inelastic scat- 
tering events: those related to charge fluctuations in 
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FIG. 1; (color online) Non-equilibrium charge susceptibility 
(full lines) and dynamical conductance G (dashed lines) 
versus applied bias, (a) Non-interaction case, with eo = 0.5 
and Pa = 0.7; (b-c) With interaction and for different trans- 
port regimes {Pa ~ 2). From off-resonant to resonant: (b) 
£o = 0.7, (c) eo ~ 0.5, (d) eo = 0.f5. Xc^ is rescaled by 
r = FaifJ-''^) = tla/Pa- Ou this scale both and G present 
the same spectral features: peaks associated with charge fluc- 
tuations in the electronic resonances. Calculations are done 
for symmetric coupling t^a = 0.15 and asymmetric potential 
drops /XL = n"'^ + eV, jj.R — fi'''^. The other parameters are 
ljJo = 0.3 and 70 = 0.21, e^ = 0. The energy parameters are 
in [eV] and G in unit of Go ~ /h. 

the electronic resonances (resonant elastic and inelastic 
transmission) and those involving off-resonant inelastic 
scattering by tunneling electrons. However the relation- 
ship derived previously for G{V) and X^^(V) is indepen- 
dent of the nature of the interaction (electron-phonon or 
electron-electron) in the central region C. 

In our model, the Hamiltonian for the region C is 

He — eod^ d + ujQa^ a + jQ^a^ + a)d^ d, (5) 

where d^ (d) creates (annihilates) an electron in the level 
eo, which is coupled to the vibration mode of energy 
llIq via the coupling constant 79. The central region 
C is connected to two (left and right) one-dimensional 
tight-binding chains via the hopping integral ioL a-nd 
ton. The corresponding lead a — L, R self-energy is 
S^(a;) = tgQ,//3Q, exp''^°("' with the dispersion relation 
w = £q + 2/3q, cos(fcct(a;)). Here, the imaginary part 
Ta — — 29mS^ is energy dependent and goes beyond 
the wideband limit, unless /3q. is much larger than any 
other paramaters. At equilibrium, the whole system 
has a well-defined unique Fermi level n'^'^. A finite bias 
V, applied across the junction, lifts the Fermi levels as 
Mi.fl + Vl.rsV. The fraction of potential drop at 

the left contact is 77^ and rjn = 77^ — 1 at the right contact 
[55], with — M-R = and r]L G [0, 1]. 



Finally the electron-phonon interaction is treated at 
the Hartree-Fock level (first order diagrams in term 
of the interaction) and is incorporated as self-energies 
sfp^'*^^"^^ (w) in the NEGF. Self-consistent calculations 
provide a partial resummation of the diagrams to all or- 
ders Bll^. 

Within this model, we calculate the dynamical con- 
ductance G{V) from Eq. p| and the NE charge suscep- 
tibility X^^(y) from Eq. M for different sets of parame- 
ters. We consider symmetric (toL — toR.) and asymmetric 
{toL 7^ ton) coupling to the leads, different strength of 
coupling to the leads, symmetric and asymmetric poten- 
tial drops at the contacts, and different transport regimes 
(off-resonant Eq ^ and resonant eo ~ m"'^)- We re- 
strict ourself here to the medium electron-phonon cou- 
pling (0.5 < 7o/wo < 1) regime which corresponds to 
realistic coupling in organic molecules. The strong cou- 
pling regime requires higher-order diagrams and more 
time consuming calculations [211 dS] . 

Figure [1] shows the NE charge susceptibility x?^(^) 
and the dynamical conductance G{V). We consider a 
symmetric coupling to the leads (toL = toR) and an 
asymmetric potential drop (77^ = 1). On this scale, both 
the conductance and the NE charge susceptibility present 
peaks at an applied bias corresponding to an electronic 
resonance: a main resonance peak close to full polaron 
shift renormalised level eo = eo — 7o /ojq , and phonon side- 
band peaks around V ^ Eq + ntuo [27 . In the NE condi- 
tions, the charge fluctuates in these electronic resonances 
whenever the bias window includes eo+nwo. Hence peaks 
are obtained in the charge susceptibility x?^(^) for these 
biases. The peaks correspond to elastic (I^ '-^ eo) and in- 
elastic {V ~ eo + nwo) resonant scattering processes. For 
the non- interacting case Fig. [ija), there is only one reso- 
nance at eo, and, as demonstrated, and G are related 
via r(/ii) beyond the wideband limit. 

Figure [2] shows that the relationship between x^^ and 
G is robust against our model parameters. It holds 
for asymmetric coupling to the leads (<ol 7^ ^oi?) and 
different fractions of potential drops at the contacts - 
see Fig 2(d). It holds for strong coupling to the leads 
~ 1^0 > 7o - see Fig 2(c); and beyond the wideband 
limit - see Fig. 2(b). It also holds when the interaction 
is modelled only with the Fock diagram - see Fig. 2(a). 
Therefore the relation between and G is not due to 
the fact that the Hartree self-energy E^^j^ is proportional 
to (ngE) (hence dv^^,, (x x?'')- 

Note that with potential drops at both contacts, 
and /ifl support a fraction of the bias, and the relation- 
ship between x^^ and G includes terms in dy Jr- How- 
ever both quantities still present the same features versus 
applied bias - see Fig. 2(d). 

On a smaller energy scale, the conductance also con- 
tains physical information for biases around excitation 
energies which goes beyond resonant transmission. In- 
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FIG. 2: (color online) Non-equilibrium charge susceptibility 
(full lines) and dynamical conductance G (dashed lines) 
versus V . Same parameters as in Fig. 1(c) except as otherwise 
stated, (a) Far beyond the wideband approximation: Pa = 
0.7. (b) Interaction at the Fock level only, (c) Strong coupling 
to the leads toa = 0.30 = loq. (A) Asymmetric coupling to the 
leads and potential drop, Iql ~ 0.07, toij — 0.15 and til = 
toR/{toL+toR) = 0.68182, and eq = 0.70. 
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FIG. 3: (color online) Derivative dvy^^ (dashed lines) and 
lETS signal dvG (full lines), dvx " 
dip feature of the lETS at V ~ cjq 
regime. Calculations done with = 1-50 ~ /3ci, eo = 0.0, 
7o = 0.195. (Bottom) Off-resonant regime for different eo- 
dvX^^ and dvG are normalised by Xc*^ and G respectively, 
(a) eo = 0.70, (b) eo = 0.99, (c) eo = 1-20. The arrows 
point the position of the electron resonance a.t V = io — uiq. 
Calculations done with the same parameters as in Fig. 1(b). 



deed, the conductance also varies at bias thresholds cor- 
responding to other inelastic scattering processes (for ex- 
ample inelastic electron tunneling). At the bias thresh- 
old y ~ Wo, the conductance increases in the ofF-resonant 
transport regime (opening of new conduction channels) 
or decreases in the resonant transport regime (electron- 
phonon backscattering) . These effects are better seen 
in the inelastic electron tunneling spectroscopy (lETS) 
as peaks or dips for the off-resonance or resonant trans- 
port regime respectively [25l [28] . The JETS is obtained 
from the second derivative of the current versus applied 
bias d^I/AV^ = AG/dV = dvG{V). In experiments, 
the lETS signal is usually given normalised by the con- 
ductance itself or by the current itself. Figure [3] shows 
the lETS signal as well as the corresponding variation 
of the non-equilibrium charge susceptibility versus ap- 
plied bias 5yx^^- O^i^ can clearly see a peak feature 
at y ~ Wo in the lETS signal, while x?^ is virtually 
featureless at the corresponding bias for both the off- 
resonant and resonant transport regimes. This means 
that these inelastic tunneling electron-phonon scattering 
processes (at ~ wq) are not related to charge fluctua- 
tions. Instead the phonon population fluctuates because 
of phonon emission induced by the tunneling electron. 
Note that the tiny features pointed by the arrows in Fig- 
ure [s] correspond to tiny peak features in both x?^ and 
G. They are related to charge fluctuations in the electron 
resonances at y = eo — wq (phonon emission by a hole). 

Discussion: We have hence shown that the non- 
equilibrium charge susceptibility and the dynamical con- 
ductance are directly related to each other, though in a 



different manner than for the equilibrium case. They 
both present features (peaks) versus the applied bias 
whenever there are charge fluctuations in the correspond- 
ing electronic resonances of the nano junction. 

Therefore we suggest that measuring both the conduc- 
tance and the NE charge susceptibility simultaneously in 
the same experiment is essential in quantum transport. It 
permits one to identify the nature of the scattering pro- 
cesses involve in the transport, i.e. processes involving 
charge fluctuation or not. This result is very important 
for the analysis of the transport properties in complex 
systems such as large single- molecule junction and does 
not involve the presence of a third gate electrode. Al- 
though our result is mostly relevant for electron-phonon 
scattering processes, it is not limited only to these pro- 
cesses. The measurement of the NE charge susceptibil- 
ity could be performed by measuring the potential drop 
around a capacitor placed in series with the nanojunc- 
tions (Vcap — e{nc} /Ccap)- One can then obtain x?^(^) 
in a similar way as the dynamical conductance G{V) is 
obtained from the current by using a lock-in set up. 
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